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I.  INTRODUCTION 


A  FORTRAN  program  was  written  to  find  the  physical  properties  of  one 

sabot  petal.  The  program  was  written  to  provide  a  simple  analytical  model  to 
determine  these  properties.  This  analytical  model  is  an  easy  alternative  to 
making  actual  physical  measurements.  Though  the  program  was  developed  for 

analyzing  a  sabot  petal,  It  can  be  used  for  any  homogeneous  partial  or  com¬ 
plete  body  of  revolution.  The  only  restriction  is  that  the  body  of  revolution 
must  be  defined  by  an  angle  0^.  as  shown  In  Figure  1. 

Section  II  provides  the  theory  for  finding  the  physical  properties  of 
one  sabot  petal.  The  program  takes  user-provided  contour  points  of  the 
sabot's  RZ  plane  profile  (Figure  1)  and  uses  them  to  divide  the  sabot  petal 
into  a  finite  number  of  partial  pieces.  The  program  then  finds  the  physical 

properties  of  each  partial  piece  and  by  using  the  correct  addition  formulas, 

adds  all  the  partial  piece  results  to  get  the  solution  for  the  total  sabot 
petal.  Other  than  the  coordinates  of  the  profile,  the  only  Information  the 
user  provides  Is  the  material  density  of  the  sabot  and  the  angle  of  the  sabot 
petal,  6j. 

In  Section  III,  the  procedure  for  transforming  the  inertia  tensor  to  any 
Cartesian  coordinate  is  described.  All  the  equations  are  provided,  but  the 
user  must  be  careful  of  two  things.  First,  the  user  must  not  confuse  skew 
lines  for  parallel  lines  In  the  translation  transformation  equations.  Second, 
the  user  must  follow  the  "right-handed  rule"  when  determining  the  sign  of  the 
rotation  angle  In  the  rotation  transformation  equations. 

Test  cases  are  found  In  Section  IV,  The  theory  Is  tested  on  a  25mm 
aeroballlstlc  model  by  comparing  the  results  found  by  physical  measurement 
machines  with  the  answers  predicted  by  the  program.  The  transformation  equa¬ 
tions  are  shown  to  work  on  a  uniform  rod. 

In  Appendix  A  there  Is  a  listing  of  the  program,  and  in  Appendix  B  an 
example  problem  is  solved  following  step-by-step  Instructions  on  the  execution 
of  the  program. 


II.  THEORY 

The  sabot  petal  Is  divided  Into  partial  pieces  as  required  to  describe 
Its  particular  geometry.  This  Is  done  by  describing  the  profile  of  the  sabot 
petal  In  the  RZ  plane  (Figure  1)  using  contour  points.  The  program  uses  the 
contour  points  to  construct  the  partial  pieces  and  then  finds  their  physical 
properties.  After  each  partial  piece's  physical  properties  are  found,  the 
physical  properties  of  the  total  sabot  petal  are  found  by  adding  the  partial 
piece's  properties  together  using  the  appropriate  addition  formulas. 

The  first  physical  property  found  Is  the  volume.  The  volume  for  any 
partial  piece  "j  '  of  the  total  sabot  petal  is:1 


v(j)  »  ///  dV  (1) 


1 


If  a  local  cylindrical  coordinate  system  is  used  (Figure  2),  Equation  (1) 
becomes: 


v(J) 


r  dedrdz 


(2) 


See  Figure  3  for  Integration  nomenclature  and  Figure  2  for  the  definition  of 
9j.  The  resulting  equation  for  the  volume  of  partial  piece  "j"  Is: 

v(j)  *  0^{m2zt3/3  +  z^b  +  b2zt }  .  (3) 


Note  that  the  upper  limit  of  r  Is  restricted  to  straight  line  segments. 
This  means  that  any  nonlinear  curve  In  the  ZR  plane  must  be  approximated  by 
straight  line  segments.  If  this  Is  necessary,  the  approximating  error  can  be 
minimized  by  making  the  number  of  approximating  segments  sufficiently  large. 
This  Is  the  only  approximation  used  In  the  program. 

The  lower  limit  of  r  In  Equation  (2)  Is  always  taken  as  zero.  The  sabot 
petal,  which  Is  a  body  with  cavities.  Is  handled  by  the  program  as  follows: 
The  outer  surfaces  of  the  sabot  define  solid  bodies  of  calculable  volume;  the 
Inner  surfaces,  which  are  the  boundaries  of  the  cavities,  define  other  solid 
bodies  whose  volumes  are  flagged  to  be  negative.  The  volume  for  the  total 
sabot  petal  Is  then  computed  as  the  sum  of  the  positive  and  negative  volumes. 
This  approach  made  the  programming  easier  and,  more  importantly,  kept  the 
required  Input  data  to  a  minimum.  Figure  4  Illustrates  the  above  procedure. 
This  procedure  Is  used  for  all  the  physical  properties  found  by  this  program. 

The  equation  for  the  volume  of  the  total  sabot  petal  Is: 

n 

V  «  L  v(j)  (4) 

M 


where  v(j)  Is  the  volume  of  partial  piece  "j" 
n  Is  the  number  of  partial  pieces. 

Because  only  a  uniform  density  Is  allowed,  the  total  mass  of  the  sabot 
petal  Is: 

M  -  Vp  (5) 


where  p  Is  the  density  of  the  material. 

With  the  total  mass  found,  the  center  of  mass  can  be  calculated.  The  x 
location  for  the  center  of  mass  of  any  body  Is  defined  as:2 


Xc  -  JUsUL 

J/Jp  dV 


(6) 


Since  the  density  Is  assumed  to  be  uniform.  It  can  be  taken  out  of  the 
Integrals  and  canceled  out.  This  leaves  a  denominator  equal  to  the  volume  of 
the  mass  in  question.  Since  the  volume  of  any  partial  piece  "j"  has  been 
determined,  the  only  Integral  that  needs  to  be  determined  to  calculate  the  x 
center  of  mass  for  any  partial  piece  "j"  Is: 


///  x  dV 


(7) 


In  a  local  cylindrical  coordinate  system.  Expression  (7)  becomes 


r2  cose  dedrdz 


(8) 


Expression  (8)  Is  evaluated  to  be: 


(2/3)s1  n6j{  (l/4)m3z^.1*  +  m2zt3b  +  (3/2)b2mzt2  +  b3zt}  (9) 


The  x  location  of  the  center  of  mass  for  partial  piece  "j”  In  the  global 
coordinate  system  Is 


(2/3)s1 n8, { (l/4)m3zt4  +  m2z  3b  +  (3/2)b2mz  2  +  b3zf} 
xc(j)  - - l - 1 - ^ - 1 - L.  (10) 


The  x  location  of  the  center  of  mass  for  the  total  sabot  petal  In  the  global 
coordinate  system  Is  the  sum  of  the  volumes  v(j)  multiplied  by  the  center  of 
mass  xc(j).  This  value  Is  then  divided  by  the  total  volume  of  the  sabot  petal 
to  give  XC,  the  x  location  of  the  center  of  mass  of  the  total  sabot  petal  In 
the  global  coordinate  system.2  This  Is  done  In  Equation  (11). 


XC 


E  v(j)xc(j) 

_ 

E 

j-1 


(11) 


The  definition  of  the  y  location  of  the  center  of  mass  is  similar  to  the 
definition  of  the  x  location  of  the  center  of  mass.  The  only  integration  that 
needs  to  be  performed  to  determine  the  y  center  of  mass  for  any  partial  piece 
•T  Is: 


3 


r2sin6  dedrdz 


(12) 


fZ^  y»!11Z+b  /'01 

I  {  l 


However, 


sine  de  *  0 


This  means  the  y  location  of  the  center  of  mass  of  all  the  "j"  partial  pieces 
Is  0.  This  result  was  expected  by  the  argument  of  symmetry  (Figure  2).  By 
replacing  all  the  "X,Ms  with  "Y*"s  In  Equation  (11)  and  using  the  result  from 
Equation  (13),  the  Y  location  of  the  center  of  mass  for  the  total  sabot  petal, 
YC,  Is  0. 

The  definition  of  the  z  location  of  the  center  of  mass  Is  the  same  as  the 
definition  of  the  x  and  y  locations  of  the  center  of  mass.  The  Integral, 
which  Is  similar  to  Expressions  (8)  and  (12),  Is: 


y»z«.  /•mz+b  /*0« 

ill, 


zr  dedrdz  . 


Expression  (14)  Is  evaluated  to  be: 

e^UM^-m2  +  (2/3)zt%b  +  (l/2)b2zt2}  (15) 

The  z  location  for  the  center  of  mass  of  partial  piece  "j"  In  the  global 
coordinate  system  Is: 


zc(j)  *  ZO  + 


ejUl/^z^m2  +  (2/3)zt3mb  +  (l/2)b2zt2} 


where  ZO  Is  the  location  of  the  local  coordinate  system's  origin  In  the  global 
coordinate  system  (Figure  2).  The  z  location  for  the  center  of  mass  for  the 
total  sabot  petal  In  the  global  coordinate  system  Is  found  the  same  way  the  x 
location  was  found. 

£  v( j )zc( j ) 

ZC  -  -  (17) 

i  ^ 

j*i 


i 


The  mass  moments  of  inertia  are  defined  as2 


Ixx  *  ///  p(y2  +  z2)  dV 


Iyy  3  ///  p(x2  +  z2)  dV 
Izz  *  ///  p(x2  +  y2)  dV 


Equations  (18),  (19)  and  (20)  for  partial  piece  "j"  in  its  local  cylindrical 
coordinate  system  are: 


1xx(j)  -  p 


rrt 


(z2r  +  r3sin28)de  dr  dz 


iyy(j)  *  p 


nmz+b  y*  8. 

L 


(z2r  +  r3cos2e)de  dr  dz 


y»zt  y-mz+b  ^-8^ 

izz(j)  *  pj  J  'J  r3  da  dr  dz 


0  ■'O  -0 


Equations  (21),  (22),  and  (23)  are  evaluated  to  be: 


txx(j) 


20  A  /m4zt5  bm3zt4  b2m2zt3 
I  \  20  4  2 

Vt5"12  .  Vt*""’  ,  8i2tib2) 

5  2  3  ) 


4  b2m2z  3  b3mzf2  b4zA 
-  + - -  + - -  + - L 


iyy(j)  *  p  (9,  + 


sin28 


A  /m4zt5  bm3zt4  b2m2zt3  b3mzt2  b4zA 
/  V  20  4  2  2  4/ 


0lZt5m2  0lzt4mb  0lZt3t)2 


5 


p9.  (m4Zf5  \ 

1zz(j)  » - | - +  m3zt4b  +  2b2m2zt3  +  2b3mzt2  +  b4zt(  .  (26) 


Equations  (24),  (25),  and  (26)  are  partial  piece  "j"'s  mass  moments  of  Inertia 
about  Its  own  local  coordinate  axes.  In  order  to  find  the  mass  moments  of 
Inertia  for  the  total  body,  all  the  partial  piece  results  must  be  found  about 
axes  through  the  centers  of  mass  of  the  partial  pieces.  This  Is  accomplished 
by  the  use  of  the  parallel  axis  theorem: 


Ixx  *  Ix'x'  +  Md2 


(27) 


where  Ix'x'  Is  the  moment  of  Inertia  about  an  axis  through  the  center  of  mass 


Ixx  Is  the  moment  of  Inertia  about  an  axis  that  is  parallel  to  an  axis 
through  the  center  of  mass 

M  Is  the  mass  of  the  object 

d  Is  the  distance  between  the  xx  and  x'x*  axes 


The  mass  moments  of  Inertia  of  partial  piece  "j"  about  axes  through  its  center 
of  mass  and  parallel  to  Its  local  system  are  derived  from  Equation  (27)  and 
the  results  from  Equations  (24),  (25),  and  (26).  They  are: 


Ix'x'(j)  -  1 xx (j )  -  pv(j){zc(j)  -  ZO}2 
ly’y'(j)  -  iyy(j)  -  pv(j){xc(j)2  +  (zc(j)  -  zo)2} 
Iz’z'(j)  =  1 zz ( j )  -  pv(j)xc(j)2  . 


(28) 

(29) 

(30) 


The  mass  moments  of  Inertia  for  the  total  sabot  petal  about  axes  through  its 
center  of  mass  and  parallel  to  the  global  coordinate  system  are 


IX'X*  «  Z  Ix'x' (j)  +  pv(j)(ZC  -  zc( j ) )2  (31) 

M 


IVY' 


iy’y'(j)  +  pv(j)[(zc  -  zc(j))2  +  (xc  -  xc(j))2] 


(32) 


(33) 


IZ'Z'  *  £  iz'z'(j)  +  pv(j)(XC  -  xc(j))2  . 

j=l 


Equations  (31),  (32),  and  (33)  were  derived  from  the  parallel  axis  theorem  for 
mass  moments  of  inertia  and  Equations  (28),  (29),  and  (30). 


The  mass  product  of  inertia  terms  are  defined  as2 


Ixy  =  ///  pXY  dV 

(34) 

Ixz  ■  J/J  pXZ  dV 

(35) 

lyz  -  III  pYZ  dv. 

(36) 

For  partial  piece  "j"  in  its  local  cylindrical  coordinate  system.  Equations 
(34),  (35)  and  (36)  reduce  to 


nmz+b  /-0J 

I  r3  cose  sine  de  dr  dz 

-6j 


(37) 


/-zt  /-mz+b  /.Oj 

1xz(j)  •pi  I  I  zr2  cose  de  dr  dz 

T)  0  re. 


(38) 


'1 


rz  t  rmz+b  /*e  i 

1yz(j)  ■  p  /  /  /  zr2  sine  de  dr  dz  . 

u  -b  -ej 


(39) 


Equations  (37),  (38)  and  (39)  are  evaluated  to  be 
ixy(j)  *  0 


fzt5m3  3  z  2b3l 

ixz(j)  *  (2/3)p  sine. -  +  -  z  4m2b  +  z  3mb2  +  - 

XL  5  4  C  C  2  J 


iyz(j )  =  o  . 


(40) 

(41) 

(42) 


Equations  (40),  (41)  and  (42)  are  the  mass  product  of  inertia  expressions  for 
partial  piece  "j"  about  its  local  coordinate  system.  Like  the  mass  moments  of 
inertia,  the  product  of  inertia  terms  need  to  be  found  about  axes  through  the 
center  of  mass  of  the  partial  piece  "j".  This  is  accomplished  by  using  the 
parallel  axis  theorem  for  products  of  inertia.2 


Ixy  =  Ix'y '  +  M  Xc  Yc  (43) 
Ixz  =  lx1 2'  +  M  Xc  Zc  (44) 
Iy2  =  ly ' 2 ’  +  M  Yc  Zc  (45) 


where  Ix'y',  Ix'2'  and  Iy ' z '  are  the  products  of  inertia  about  the  center  of 
mass 

Ixy,  1x2  and  Iy2  are  the  products  of  inertia  about  a  set  of  axes 
parallel  to  the  center  of  mass 

Xc,  Yc  and  Zc  are  the  x,  y  and  2  locations  of  the  center  of  mass  in  an 
XYZ  coordinate  system. 


Equations  (44)  and  (41)  are  used  to  generate  the  equation  for  Ixz  for  partial 
piece  "j"  about  its  own  center  of  mass: 


1x'2'(j)  =  ixz(j)  -  pv(j)xc(j)(2c(j)  -  ZO)  .  (46) 


1x2  for  the  total  sabot  petal  about  axes  through  its  center  of  mass  and 
parallel  to  the  global  coordinate  system,  is 


n 

I X ' Z '  -  £  i x 1 z 1 ( j )  +  pv(j)(xc(j)  -  XC)(zc(j)  -  ZC)  .  (47) 

j=l 


Equation  (47)  was  arrived  at  by  the  use  of  the  parallel  axis  theorem  for  prod¬ 
ucts  of  inertia  and  Equation  (46).  Because  the  y  center  of  mass  Is  zero  for 
every  "j"  partial  piece  (Equation  (13))  and  1xy(j)  is  zero  for  every  "j"  par¬ 
tial  piece  (Equation  (40)),  by  Equation  (43),  ix'y'(j)  for  every  "j"  partial 
piece  is  zero.  The  total  sabot  petal's  Ix'y',  IX'Y',  is  zero  by  the  same 
argument.  The  same  argument  holds  for  Iy ' z '  for  the  total  sabot  petal;  there¬ 
fore  IY'Z'  is  zero. 


III.  INERTIA  TENSOR  TRANSFORMATIONS 


The  program  calculates  the  Inertia  tensor  about  a  set  of  axes  that  go 
through  the  center  of  mass  of  the  sabot  petal.  This  section  explains  how  to 


transform  the  inertia  tensor  properly  from  this  coordinate  system  to  any  other 
Cartesian  coordinate  system.  There  are  two  transformation  procedures:  trans¬ 
lation  and  rotation.  Either  procedure  can  be  done  Independently  or  they  can 
be  done  in  successive  steps.  The  equations  of  this  section  make  use  of  the 
results  found  by  the  program  as  much  as  possible  for  ease  of  use.  Care  should 
be  taken  If  these  equations  are  to  be  used  for  anything  but  this  purpose. 

Translation  of  the  Inertia  tensor  requires  the  use  of  both  types  of 
parallel  axis  theorems  used  in  Section  II  of  this  report.  Remember  to  use 
consistent  units  and  make  sure  the  axes  are  parallel.  The  mass  moment  of 
Inertia  terms  are  translated  by  Equation  (27).  With  the  results  found  by  the 
program,  the  general  Equation  (27)  Is  rewritten  as  the  following  equations. 


I'XX  *  IX'X'  +  M[(ZC  -  Z'O)2  +  (-Y'O)2]  (48) 


I'YY  *  IY'Y'  +  M[(XC  -  X'O)2  +  (ZC  -  Z'O)2]  (49) 


I'ZZ  *  IZ'Z*  +  M[(XC  -  X'O)2  +  (-Y'O)2]  (50) 


where  X'O,  Y'O  and  Z'O  are  the  X,  Y  and  Z  locations  of  the  desired  coordinate 
system's  origin  In  the  global  coordinate  system. 

For  the  product  of  Inertia  terms.  Equations  (43),  (44)  and  (45)  are 
rewritten  as: 


i'xy  =*  m(xc  -  x'o)(-Y'n) 

(•51) 

rxz  =  ix' z'  +  m(xc  -  x'o)(zc  -  z*o) 

(52) 

I'YZ  =  M(-Y'O) (ZC  -  Z'O)  . 

(53) 

Note  that  all  the  values  in  Equations  (48-53)  are  found  by  the  program  except 
X'O,  Y'O  and  Z'O,  which  are  user-defined. 

Only  one  translation  is  required  to  translate  the  coordinate  system  of 

the  Inertia  tensor  found  by  the  program  to  any  other  coordinate  system  in 

space.  This  Is  Important  for  the  user  to  know  because  Equations  (48-53)  are 
only  valid  for  the  case  where  the  moments  and  products  of  inertia  are  known 

about  axes  that  go  through  the  center  of  mass  of  the  body.  This  means  the 

user  cannot  take  the  results  from  one  translation  and  use  Equations  (48-53)  to 
make  another  translation.  In  general,  only  one  translation  is  allowed. 

Rotation  of  the  coordinate  system  Is  allowed  and  one  proper  way  to  do 
this  is  presented.  The  equation  to  rotate  the  inertia  tensor  is:3 
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CT']  =  [B][T][B] 


(54) 
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where: 


tT']  is  the  Inertia  tensor  about  the  new  coordinate  system 
T]  Is  the  Inertia  tensor  about  the  old  coordinate  system 


[T]  Is  the  Inertia  tensor  about  the  old  coordinate  system 
[B]  is  the  transformation  rotation  tensor 

[B]‘l  Is  the  inverse  of  the  transformation  rotation  tensor. 


Equation  (54)  allows  rotation  about  a  point,  but  for  simplicity  only  the  rota¬ 
tion  about  any  of  the  coordinate  axes  will  be  found.  This  may  require  more 
than  one  rotation  to  get  to  the  desired  coordinate  system,  but  this  Is  allowed 
as  long  as  the  procedure  Is  followed  as  directed.  Before  Equation  (54)  can  be 
solved,  the  B  tensor  must  be  determined.  Figure  5  and  Equations  (55-63)  show 
how  to  transform  the  coordinates  of  an  EFG  coordinate  system  to  coordinates  In 
an  E'F'G'  coordinate  system.  The  transformation  is  a  rotation  about  the  G 
axis;  therefore  G'  *  G. 


d^  =  ftane 


4  *  e  +  dj 


e'  ■  Ljcose 


e*  *  ecosfl  +  fslne 


d2  *  etane 


L2  *  f  -  d2 


f'  =  L2C0S9 


f'  ■  -eslne  +  fcose 


g'  a  g  . 


In  matrix  form,  the  rotation  transformation  from  the  EFG  coordinate 
system  to  the  E'F'G'  coordinate  system  is: 


cose  sine 


-sine  cose 


mmvsiffffig 


nmmnn 


Therefore 


[8] 


cose  sine  0‘ 

-sine  cose  0 

LOO  1, 


(65) 


A  property  of  a  transformation  tensor  that  rotates  one  Cartesian  coordinate 
system  to  another  is:4 


[B]-1  *  [B]T 


(66) 


where  the  "T“  superscript  denotes  a  transposed  matrix. 

Equations  (65)  and  (66)  allow  Equation  (54)  to  be  written  In  matrix  form 


as: 


“I  EE ' 

IEF' 

IEG1' 

*  COS0 

si  n0 
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where: 

0  Is  the  angle 

of  rotation 

shown 

In 

Figure 

5 

• 

(67) 


The  equations  that  result  from  the  the  solution  of  Equation  (67)  are 

I  EE'  -  IEEcos26  +  IFFs1n28  +  IEFs1n20  (68) 

IFF'  =  IEEs1n20  +  IFFcos2e  -  IEFs1n20  (69) 

IGG'  *  IGG  (70) 

IEF*  «  IFE'  -  .5( IFF  -  IEE)sin26  +  IEF(cos26  -  s1n2e)  (71) 

I  EG*  -  IGE'  »  IEGcose  +  IFGsIne  (72) 

IFG *  -  IGF'  »  IFGcose  -  IEGsIne  .  (73) 

Equations  (68-73)  are  generalized  equations  and  to  use  them  In  the  XYZ  coordi¬ 
nate  system  make  the  following  substitutions.  For  a  rotation  about  the 


X  axis 

Let  Y  *  E 
Z  «  F 
X  «  6 


Y  axis 

Let  Z  »  E 
X  *  F 
Y  *  G 
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Z  axis 

Let  X  =  E 
Y  *  F 
Z  =  G 


In  Equations  (68-73). 

The  above  Instructions  must  be  followed  to  Insure  a  proper  transformation 
from  one  right-handed  Cartesian  coordinate  system  to  another.  As  many  rota¬ 
tions  as  required  are  allowed  using  Equations  (68-73).  If  a  translation  and 
rotation  of  the  inertia  tensor  are  to  be  performed,  translate  by  Equations 
(48-53),  then  rotate  as  many  times  as  necessary  by  Equations  (68-73).  If  a 
rotation  of  the  axis  system  is  performed.  Equations  (51)  and  (53)  are  In 
general  no  longer  valid. 


IV.  TEST  CASES 

The  purpose  of  the  program  was  to  provide  an  alternative  to  making  actual 
physical  measurements  to  determine  the  physical  properties  of  one  sabot  petal. 
The  test  of  the  program  was  to  compare  the  physically  measured  properties  with 
properties  found  by  the  program.  The  sabot  petal  used  for  the  test  is  from  a 
25mm  aeroballistic  model.  The  petal  Is  made  from  aluminum  and  Is  one  quarter 
of  a  complete  sabot  package.  It  is  shown  in  Figure  6. 

The  sabot  petal  mass,  center  of  mass  and  mass  moments  of  inertia  were 
found  by  physical  measurement  machines  located  at  the  Transonic  Range  facility 
of  Aberdeen  Proving  Ground.  These  machines  use  precision  air  bearings  and 
have  an  accuracy  of  0.03%.  The  program  determined  the  properties  using  22 
contour  points.  The  results  for  this  particular  case  have  a  maximum  differ¬ 
ence  of  4.8%  when  compared  with  the  results  from  physical  measurement  equip¬ 
ment.  if  203  points  are  used  to  more  accurately  describe  the  buttress  grooves 
and  circular  arcs,  the  maximum  difference  drops  to  2.9%.  The  comparison  of 
the  three  cases  Is  shown  in  Figure  6.  This  shows  a  trend  towards  smaller 
percentage  differences  when  the  number  of  contour  points  Is  increased.  If  the 
theory  is  correct,  this  trend  should  occur. 

An  example  of  the  use  of  the  transformation  equations  is  provided.  A 
uniform  rod  Is  divided  Into  six  pieces  as  shown  In  Figure  7.  The  inertia 
properties  of  each  piece  are  determined  by  the  program.  If  the  transformation 
equations  are  correct,  the  total  rod's  Inertia  tensor  about  its  center  of  mass 
can  be  found  by  transforming  all  the  partial  piece  results  to  the  location  of 
the  center  of  mass  of  the  total  rod  and  adding  them  together. 

TABLE  1  shows  the  physical  properties  calculated  by  the  program  for  each 
partial  piece.  These  results  are  then  translated  using  Equations  (48-53). 
The  translated  results  are  shown  In  TABLE  2.  The  translated  results  are  then 
rotated  using  Equations  (68-73)  and  the  results  are  shown  In  TABLE  3.  With 
these  transformations  completed,  all  the  partial  pieces'  Inertia  tensors  are 
known  about  the  same  set  of  axes,  a  set  whose  origin  Is  the  center  of  mass  of 
the  total  rod.  This  allows  the  results  found  In  TABLE  3  to  be  added  and  the 
sums  should  equal  the  Inertia  terms  for  the  total  rod. 


TABLE  3.  Tensor  Rotations. 


1 

2 

3 

4 

5 

6 

Sum 

1-6 

Ixx 

5.257 

.890 

3.384 

2.332 

9.792 

5.576 

27.231 

iyy 

5.257 

.868 

3.516 

2.267 

9.704 

5.621 

27.233 

izz 

.158 

.026 

.238 

.159 

.740 

.264 

1.585 

Ixy 

0.000 

.007 

.038 

-.057 

.026 

-.013 

.001 

Ixz 

-.538 

.036 

.552 

.202 

.502 

-.754 

0.000 

Iyz 

0.000 

-.134 

-.148 

.349 

.135 

-.202 

0.000 

The  Inertia  terms  for  a  homogeneous  right  circular 
the  following  equations:2 

cylinder  are  found  by 

Ixx  -  Iyy  -  (l/12)M(3rz  +  h2) 

(74) 

Izz  »  (l/2)Mr2 

(75) 

Ixy  ■  Ixz  ■  Iyz  «  0 

(76) 

where  M  Is  the  mass  of  the  rod,  r  Is  the  radius  of  the  circular  cross  section 
and  h  Is  the  height  of  the  cylinder. 

The  Inertia  terms  for  the  total  rod  are 


Ix'x' 

-  Iy'y' 

-  27.235  lb-1n2 

(77) 

Iz'z' 

«  1.587 

1  b-1n2 

(78) 

Ix'y* 

«  Ix'z' 

-  ly'z'  -  0 

(79) 

If  the  transformation  equations  are  correct,  the  expected  results  found  by 
Equations  (77-79)  should  equal  the  results  found  In  the  sum  row  In  TABLE  3. 
It  has  been  shown  that  for  this  particular  problem,  the  transformation  equa¬ 
tions  work.  The  small  differences  in  the  expected  results  are  due  to  round¬ 
off  error. 


V.  CONCLUSIONS 

Inertia  properties  are  sometimes  Important  to  know.  Shapes  like  sabot 
petals  usually  make  It  difficult  to  calculate  Inertia  terms  analytically.  The 
alternative  to  an  analytical  calculation  Is  to  make  actual  physical  measure¬ 
ments.  Access  to  physical  measurement  equipment  Is  sometimes  difficult  and 
getting  the  necessary  Information  can  take  some  time.  The  Intention  of  this 
program  Is  to  provide  this  Information  quickly  and  accurately. 

The  Inertia  properties  of  an  object  are  defined  about  coordinate  axes. 
Because  the  inertia  Information  may  be  needed  about  axes  other  than  the 
defined  axes,  transformation  equations  were  derived.  The  report  Includes  all 
the  equations  needed  to  translate  and  rotate  the  Inertia  tensor  to  the 
required  coordinate  system.  The  equations  were  shown  to  work  In  a  selected 
example.  The  equations  should  make  the  process  of  coordinate  transformation 
easier  for  the  user. 
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Figure  3.  Integration  nomenclature. 
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Figure  4.  Lower  limit  of  r  explanation. 
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Figure  6.  Machine  vs  program  test  case. 
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Figure  7.  Inertia  tensor  transformation  example 
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I FG ’  mass  product  of  inertia  about  the  F'  and  G'  axes 

IGG  mass  moment  of  inertia  about  the  G  axis 

IGG ’  mass  moment  of  inertia  about  the  G'  axis 

ixx(j)  partial  piece  "j"'s  mass  moment  of  inertia  about  its  own  local  x 
axis 

Ix'x'(j)  partial  piece  "j"‘s  mass  moment  of  Inertia  about  an  axis  through  its 
center  of  mass  and  parallel  to  its  own  local  x  axis 

Ixx  a  body's  mass  moment  of  inertia  about  the  X  axis  in  an  XYZ 

coordinate  system 

Ix'x'  a  body's  mass  moment  of  inertia  about  an  axis  through  its  center  of 
mass  and  parallel  to  the  X  axis  in  an  XYZ  coordinate  system 

IX'X'  the  total  sabot's  mass  moment  of  inertia  about  an  axis  through  its 
center  of  mass  and  parallel  to  the  X  global  axis 

I 'XX  mass  moment  of  inertia  of  the  total  sabot  about  a  user-defined  X 
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ixy(j)  partial  piece  "j'"s  mass  pi  oduct  of  inertia  about  its  own  local  x 
and  y  axes 

Ix'y'(j)  partial  piece  "j'"s  mass  product  of  inertia  about  axes  through  its 
center  of  mass  and  parallel  to  its  own  local  x  and  y  axes 

Ixy  a  body's  mass  product  of  inertia  about  the  X  and  Y  axes  in  an  XYZ 

coordinate  system 


Ix'y'  a  body's  mass  product  of  Inertia  about  axes  through  its  center  of 
mass  and  parallel  to  the  X  and  Y  axes  In  an  XYZ  coordinate  system 

IX'Y'  the  sabot's  total  mass  product  of  inertia  about  axes  through  its 
center  of  mass  and  parallel  to  the  global  X  and  Y  axes 

I'XY  mass  product  of  Inertia  of  the  total  sabot  about  user-defined  X  and 
Y  axes 


1xz(j)  partial  piece  "j'"s  mass  product  of  inertia  about  its  own  local  x 
and  z  axes 

ix'z'(j)  partial  piece  "j‘"s  mass  product  of  inertia  about  axes  through  its 
center  of  mass  and  parallel  to  its  own  local  x  and  z  axes 


a  body's  mass  product  of  inertia  about  the  X  and  Z  axes  in  an  XYZ 
coordinate  system 
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I x ' 2 '  a  body's  mass  product  of  Inertia  about  axes  through  its  center  of 
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center  of  mass  and  parallel  to  the  global  X  and  Z  axes 

I'XZ  mass  product  of  inertia  of  the  total  sabot  about  user  defined  X  and 
Z  axes 

1yy(j)  partial  piece  "j'"s  mass  moment  of  Inertia  about  Its  own  local  y 
axis 

ly'y'U)  partial  piece  "j'"s  mass  moment  of  Inertia  about  an  axis  through  Its 
center  of  mass  and  parallel  to  Its  own  local  y  axis 

Iyy  a  body's  mass  moment  of  Inertia  about  the  Y  axis  In  an  XYZ 

coordinate  system 

Iy'y'  a  body's  mass  moment  of  Inertia  about  an  axis  through  Its  center  of 
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center  of  mass  and  parallel  to  the  Y  global  axis 

I ' YY  mass  moment  of  Inertia  of  the  total  sabot  about  a  user-defined  Y 
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center  of  mass  and  parallel  to  Its  own  local  y  and  z  axes 

lyz  a  body's  mass  product  of  inertia  about  the  Y  and  Z  axes  in  an  XYZ 

coordinate  system 

ly'z'  a  body's  mass  product  of  Inertia  about  axes  through  its  center  of 
mass  and  parallel  to  the  Y  and  Z  axes  In  an  XYZ  coordinate  system 

I Y ' Z '  the  sabot's  total  mass  product  of  inertia  about  axes  through  Its 

center  of  mass  and  parallel  to  the  global  Y  and  Z  axes 

I'YZ  mass  product  of  Inertia  of  the  total  sabot  about  user-defined  Y  and 
Z  axes 

izz(j)  partial  piece  "j"'s  mass  moment  of  inertia  about  its  own  local  z 
axi  s 

Iz'z'(j)  partial  piece  "j'"s  mass  moment  of  inertia  about  an  axis  through  its 
center  of  mass  and  parallel  to  its  own  local  z  axis 
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a  body's  mass  moment  of  inertia  about  the  Z  axis  in  an  X YZ 
coordinate  system 

a  body's  mass  moment  of  inertia  about  an  axis  through  its  center  of 
mass  and  parallel  to  the  Z  axis  in  an  XYZ  coordinate  system 

the  total  sabot's  mass  moment  of  inertia  about  an  axis  through  its 
center  of  mass  and  parallel  to  the  Z  global  axis 

mass  moment  of  inertia  of  the  total  sabot  about  a  user  defined  Z 
axis 

used  as  an  index  to  denote  any  partial  piece  used  to  describe  the 
total  body 

local  variable;  dr/dz;  used  with  b  to  define  the  upper  bound  of 
integration  on  r 

mass 

number  of  partial  pieces 
local  cylindrical  coordinate 
Inertia  tensor 

Inertia  tensor  after  a  coordinate  transformation 

partial  piece  "j"'s  volume 

the  sabot's  total  volume 

local  Cartesian  coordinate 

Cartesian  coordinate 

X  location  of  a  user  defined  system's  origin  in  the  global 
coordinate  system 

partial  piece  "j'"s  X  value  of  the  center  of  mass  in  the  global 
coordinate  system 

a  body's  X  value  of  center  of  mass 

the  total  sabot's  X  value  of  the  center  of  mass  in  the  global 
coordinate  system 

local  Cartesian  coordinate 


Y 


Cartesian  coordinate 
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Y'O 

yc(j) 

Yc 

YC 


z 

z 


zc(j) 


Y  location  of  a  user  defined  system's  origin  In  the  global 
coordinate  system 

partial  piece  "j"'s  Y  value  of  the  center  of  mass  In  the  global 
coordinate  system 

a  body's  Y  value  of  center  of  mass 

the  total  sabot's  Y  value  of  the  center  of  mass  in  the  global 
coordinate  system 

local  Cartesian  and  cylindrical  coordinate 

Cartesian  and  cycllndrical  coordinate 

partial  piece  "j'"s  Z  value  of  the  center  of  mass  In  the  global 
coordinate  system 


ZO  Z  location  of  the  partial  piece  "j'"s  origin  In  the  global 

coordinate  system 


Z'O  Z  location  of  a  user  defined  system's  origin  in  the  global 

coordinate  system 

Zc  a  body's  Z  value  of  center  of  mass 

ZC  the  total  sabot's  Z  value  of  the  center  of  mass  In  the  global 

coordinate  system 

zt  local  variable;  z  length  of  partial  piece  "j"  and  limit  of 

1  Integration  on  z 


Greek  Symbols 

de  angular  differential 

e  cylindrical  angular  coordinate 

9j  limit  of  Integration  in  local  coordinate  system;  see  Figure  2 

0T  total  angular  measurement  on  one  sabot  petal;  see  Figure  1;  note 

that  20j  »  0j 

p  density 

l  summation  sign 


APPENDIX  A.  LISTING  OF  THE  PROGRAM 


General  Notes 

1.  Program  is  written  in  FORTRAN  77 

2.  Arrays  Z  and  R  hold  the  body  contour  points.  The  dimension  of  the 
arrays  Is  set  at  100.  This  allows  for  100  contours  to  describe  the  body. 
This  can  be  Increased  if  necessary. 

3.  Array  D(j,k)  holds  all  the  information  about  each  partial  piece. 

j=number  of  the  partial  piece 

k*l  holds  the  volume  of  partial  piece  "j" 

k*2  holds  the  XC  in  the  global  system  of  partial  piece  “j" 

k=3  holds  the  ZC  in  the  global  system  of  partial  piece  "j" 

k=4  holds  the  IX* X*  of  partial  piece  "j" 

k*5  holds  the  IY'Y*  of  partial  piece  "j" 

k=6  holds  the  IZ'Z*  of  partial  piece  "j" 

k=7  holds  the  IX ’ Z '  of  partial  piece  "j". 

The  dimension  of  the  D  array  is  set  at  (100,7).  If  more  than  100  partial 
pieces  are  needed  to  describe  the  body,  increase  “j"  only.  Leave  "k"  set  at 

7. 

4. ~  The  computer  for  which  the  FORTRAN  program  was  written  uses  I/O  unit 
"5"  to  read  information  from  the  screen  and  "6"  to  write  information  to  the 
screen.  If  this  Is  not  the  case  on  the  system  where  this  program  is  going  to 
be  placed,  the  correct  read  and  write  numbers  will  have  to  be  exchanged  for 
"5"  and  "6"  in  the  program  listing. 

5.  The  Important  variable  names  in  the  program  are  described. 

All  "TOT"  variables  like  "VOLTGT"  and  "IXZTOT"  are  the  values  for  the 
total  body. 

All  "PAR"  variables  like  "VOLPAR"  and  "XCPAR"  are  partial  values  for 
the  total  body. 

NUMREG  *  number  of  partial  pieces. 

NUWTS  *  number  of  contour  points  to  describe  the  body. 

THDEG  *  0j  in  degrees.  See  Figure  1. 

TH  *  6^  in  radians.  See  Figure  2. 

RHO  =  density  of  the  material.  Units  are  defined  by  user's  entries. 

B  =  r  Intercept  of  the  top  boundary  of  a  partial  piece  in  its  own 
coordinate  system.  Units  are  defined  by  user's  entries.  It  is  b  in  Figure  3. 


M  *  slope  (dr/dz)  of  the  top  boundary  of  a  partial  piece  in  its  own 
coordinate  system.  It  is  m  in  Figure  3. 

Z  *  length  of  a  partial  piece.  Units  are  defined  by  user's  entries. 
It  Is  zt  In  Figure  3. 

"T"  variables  like  T1  or  T5  are  used  for  clarity.  They  represent 
terms  in  the  large  physical  property  equations. 


*  PROGRAM  TO  COMPUTE  THE  PHYSICAL  PROPERTIES  OF  ONE  SABOT  PETAL 

* 

*  DIMENSION  AND  DATA  ENTRY 

★ 

DIMENSION  0(100,7) ,R(100) ,Z(100) 

REAL  IXXTOT.IYYTOT.IZZTOT.IXZTOT 
CHARACTER  FNAMEM5 
PI-3.1415927 

WRITE(6 ,*) 1  ENTER  FILE  NAME' 

READ(5,101)FNAM£ 

0PEN(10,FILE=FNAME, STATUS- 'OLD') 

1-1 

4  READ(10,*)Z(I),R(I) 

IF(I.EQ.1)G0T08 

IF((Z(I).EQ.Z(1)).AN0.(R(I).EQ.R(1)))G0T09 

8  I--I+1 
G0T04 

9  READ(10,*)RH0,THDEG 
CLOSE (10) 

TH«THDEG*P 1/360. 

3  NUMPTS-I-1 

NUMREG-0 

* 

*  FINOS  THE  PROPERTIES  OF  ALL  THE  PARTIAL  PIECES  AND  STORES  THEM  IN 

*  THE  "D"  ARRAY 

* 

DO  5  K-l.NUMPTS 
IF(Z(K).EQ.Z(K+1))G0T05 
NUMREG-NUMREG+1 
L-NUMREG 

CALL  V0LCEN(R(K+1) ,R(K) ,Z(K+1) ,Z(K) ,D(L,1) ,D(L,2) ,D(L,3) ,TH) 

IF  (0(L,1).EQ.0.)G0T05 

CALL  RMXANY(R(K+1),R(K),Z(K+1),Z(K),D(L,1),D(L,2),D(L,3), 

:D(L,4) ,n(L,5) ,TH,RH0) 

CALL  RM0IZZ(R(K+1) ,R(K) ,Z(K+1) ,Z(K) ,D(L,1) ,n(L,2) ,D(L,6) .TH.RHO) 

CALL  RM0IXZ(R(K+l)fR(K),Z(K+l),Z(K),n(L,l),D(L,2),D(L.3), 

:D(L, 7), TH.RHO) 

5  CONTINUE 
* 

#  FINOS  THE  VOLUME, X  AMD  Z  LOCATION  OF  THE  CENTER  OF  MASS  FOR  THE  SABOT  PETAL 

* 

VOLPAR-O. 

XCPAR-O. 


ZCPAR-O. 

DO  6  I-l.NUMREG 

VOLPAR*VOLPAR+D(I,l) 

XCPAR»XCPAR+D(1,1)*D(I,2) 

6  ZCPAR»ZCPAR+D( I ,1)*D( 1,3) 

*  (EQ.4) 

VOLTOT-VOLPAR 

*  (EQ.ll) 

XCTOT-XCPAR/VOLTOT 

*  (EQ.17) 

ZCTOT-ZCPAR/VOLTOT 

* 

*  CALCULATES  IX*  X*  ,1 Y' Y* ,IZ*  Z* ,  AND  IX*  Z*  FOR  THE  SABOT  PETAL 

* 

DO  7  I-l.NUMREG 

*  (EQ.31) 

IXXTOT*IXXTOT+D( I ,4)+(ZCT0T-D{ I ,3) )*(ZCTQT-D( 1,3) )*RHO*D( 1,1) 

*  (EQ.32) 

I YYTOT*I YYTOT+D( 1 ,5)+( ( XCTOT-D ( I ,2) )* (XCTOT-D ( I ,2) )+ 

( ZCTOT-D ( I,3))*( ZCTOT-D (1,3))) *RHO*D (1,1) 

*  (EQ.33) 

IZZTOT»IZZTOT+D( I ,6 ) + ( XCTOT-D ( I ,2) )*(XCTOT-D( I ,2) )*RHO*D( 1,1) 

*  (EQ.47) 

7  IXZTOT=IXZTOT+D( I ,7)+(D( I ,2)-XCTOT)*(D( 1 ,3)-ZCTOT)*RHO*D( 1,1) 

* 

*  DATA  OUTPUT  AND  FORMAT  STATEMENTS 

* 

WR1TE(6,102)FNAME 
WRITE(6,103)RH0 
WRITE(6,104)THDEG 
WRITE(6 ,105)V0LT0T*RH0 
WRITE(6 ,106)XCT0T 
WRITE(6,107) 

WRITE(6,108)ZCT0T 
WRITE(6,109) IXXTOT 
WRITE(6 ,110) I YYTOT 
WRITE (6 ,11 1) IZZTOT 
WRITE(6 ,112) 

WRITE (6,11 3) IXZTOT 
WRITE(6,114) 

101  FORMAT (A4 5) 

1D2  FORMAT ( '  '.'FILE  NAME:  ' ,A45) 

103  FORMAT('  '.'DENSITY:  '  .G11.4E2) 

104  FORMAT! '  ',  'THETA  T:  \F6.2,'  DEGREES',/) 

105  FORMAT( '  '.'MASS:  '.G12.6E2,/) 

106  FORMAT! '  ','XC:  ' .G12.6E2) 

107  FORMAT!'  ','YC:  0.') 

108  FORMAT( '  ','ZC:  ' .G12.6E2,/) 

109  FORMAT( '  ' ,  ‘  IX"  X"  :  '  .G12.6E2) 

110  FORMAT ( '  ','IY"Y":  '  .G12.6E2) 

111  FQRMAT( '  ','IZ"Z":  '  .G12.6E2,/) 

112  FORMAT ( 1  ' , ' I X  "  Y  "  :  0.') 

113  FORMAT  '  '  , 1 IX"  Z"  :  '  .G12.6E2) 

114  FORMAT ( '  ' , ' I Y'  *2*  ' :  0.') 


END 
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VOLCEN  CALCULATES  THE  VOLUME  AND  CENTER  OF  MASS  FOR  ANY  PARTIAL  PIECE 

SUBROUTINE  VOLCEN(R2 ,R1 , Z2 , Z1 ,VOL ,XC ,ZC  ,TH) 

REAL  M 
B«RI 

M*(R2-R1) /(Z2-Z1) 

Z-Z2-Z1 

Tl»M*M*Z*Z*Z/3. 

T2-B*M*Z*Z 

T3*B*B*Z 

(EQ.3) 

VOL*TH*(Tl+T2+T3) 

IF(VOL.EQ.O)RETURN 

T4«M*M*M*Z*Z*Z*Z/4. 

T5»B*M*M*Z*Z*Z 

T6-B*B*Z*Z*M*3./2. 

T7«B*B*B*Z 

(EQ.IO) 

XC»2./3.*$IN(TH)*(T4+T5+T6+T7)/VOL 

T8-M*M*Z*Z*Z*Z/4. 

T9-B*M*Z*Z*Z*2./3. 

T10«8*B*Z*Z/2. 

(EQ.16) 

ZC»Z1+(TH*(T8+T9+T10) ) /VOL 

RETURN 

END 

RMXANY  CALCULATES  IX' X'  AND  IY'Y'  FOR  ANY  PARTIAL  PIECE 

SUBROUT I NE  RMX AN Y ( R2 , R1 , ZZ , Z1 , VOL , XC , ZC , I X X , I Y Y , TH , RHO ) 

REAL  M.IXX.IYY 
B*R1 

M»(R2-R1)/(Z2-Z1) 

Z-Z2-Z1 

Tl-M*M*M*M*Z*Z*Z*Z*Z/20. 

T2«B*M*M*M*Z*Z*Z*Z/4. 

T3-B*B*M*M*Z*Z*Z/2. 

T4»B*B*B*M*Z*Z/2. 

T5-B*B*8*B*Z/4. 

T6«Z*Z*Z*Z*Z*M*M/5. 

T7-Z*Z*Z*Z*M*8/2. 

T8-Z*Z*Z*B*B/3. 

X»TH-.5*SIN(2*TH) 

(EQ.24) 

IXX-RHO*(X*(Tl*T2+T3+T4+T5)+TH*(T6+T7+T8)) 

(EQ.28) 

IXX»IXX-(ZC-Z1)*(ZC-Z1)*V0L*RH0 

X«TH+.5*SIN(2*TH) 

(EQ.25) 

I YY»RHO*(X*(Tl+T2+T3+T4+T5)+TH#(T6+T7+T8) ) 

(EQ.29) 

I YY«IYY-(XC*XC+(ZC-Z1)*(ZC-Z1) )*VOL*RHO 
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RETURN 

END 


*  RMOIZZ  CALCULATES  IZ'Z*  FOR  ANY  PARTIAL  PIECE 

* 

SUBROUTINE  RM0IZZ(R2,R1,Z2,Z1,V0L,XC,IZZ,TH,RH0) 

REAL  M.IZZ 

B-Rl 

M*«(R2-R1)  /(Z2-Z1) 

Z-Z2-Z1 

Tl«M*M*M*M*Z*Z*Z*Z*Z/5 . 

T2»8*M*M*M*Z*Z*Z*Z 

T3«2.*B*B*M*M*Z*Z*Z 

T4-B*B*B*M*Z*Z*2. 

T5-B*B*B*B*Z 

*  (EQ.26) 

IZZ«RH0*TH*.  5*(T1+T2  •■'•T3+T4+T5 ) 

*  (EQ.30) 

IZZ»IZZ-VOL*RHO*XC*XC 

RETURN 

END 

* 

*  CALCULATES  IX'Z*  FOR  ANY  PARTIAL  PIECE 

* 

SUBROUT I NE  RMO I XZ ( R2 . R1 , Z2 , Z1 , VOL , XC , ZC , I XZ ,TH , RHO ) 

REAL  M.IXZ 

B«R1 

M-(R2-R1)/(Z2-Z1) 

Z-Z2-Z1 

Tl«Z*Z*Z*Z*Z*M*M*M/5. 

T2»Z*Z*Z*Z*M*M*B*.75 

T3-Z*Z*Z*M*B*B 

T4-Z*Z*B*B*B/2. 

*  (EQ.41) 

IXZ-RHO*2./3.*SIN(TH)*{Tl+T2+T3+T4) 

*  (EQ.46) 

I XZ- IXZ-RHO*VOL*XC* ( ZC- Z1 ) 

RETURN 

END 
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APPENDIX  B.  EXAMPLE  PROBLEM 


The  physical  properties  of  a  sabot  with  the  a  ZR  profile  shown  In  Figure 
8-1  are  to  be  determined.  The  sabot  petal  is  made  of  aluminum  and  Is  one- 
third  of  a  complete  sabot  package. 


STEP  1  FIND  THE  COORDINATES  OF  THE  CONTOUR  POINTS  IN  THE  ZR  PLANE. 

SEE  FIGURE  B-l. 

STEP  2  CREATE  A  DATA  FILE  OF  THE  FORM  OF  FIGURE  B-2. 

This  file  contains  the  coordinates  found  In  STEP  1.  The  data  file 
entries  must  follow  the  contour  of  the  body  in  a  clockwise  manner  until 
the  first  point  is  reached.  The  first  point  is  then  entered  again.  This 
Is  necessary  because  this  signals  to  the  computer  that  the  profile  is 
completely  entered.  The  choice  of  the  first  point  is  arbitrary. 

The  last  record  of  the  data  file  has  two  values.  The  first  value  is 
the  density  and  the  second  value  is  the  angle  9y.  See  Figure  1.  For  the 

example  problem,  6-j-  -120.  Aluminum,  the  material  of  the  example  problem, 

has  a  density  of  .00314  slug/inch.3  Please  note  that  the  length  units  of 
density  are  the  same  as  the  length  units  used  to  describe  the  ZR  profile. 

STEP  3  EXECUTE  THE  PROGRAM 

The  only  entry  required  is  the  name  of  the  file  created  in  STEP  2. 
STEP  4  SOLUTION 

The  solution  Is  written  to  the  screen  of  the  terminal.  If  a  terminal 
Is  not  used,  the  program  will  have  to  be  modified  by  the  user.  The 
example  problem's  solution  Is  shown  In  Figure  B-3.  The  units  are  con¬ 
sistent  with  the  units  entered  In  STEP  2.  For  the  example  problem  the 
mass  Is  In  slugs,  the  X  and  Z  locations  of  the  center  of  mass  are  In 
Inches  from  the  origin  of  the  global  coordinate  system  (Figure  1)  and  the 
moments  and  product  of  inertia  are  In  slugs  inches.2  Remember  the 
moments  of  inertia  are  determined  about  axes  through  the  center  of  mass 
and  parallel  to  the  global  coordinate  system. 
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FILE  NAME:  EPOF 
DENSITY:  0.3140E-02 
THETA  T:  120.00  DEGREES 

•MASS:  0.749573E-02 

XC:  0.826570 
YC:  0. 

ZC:  1.71011 


IX'X'  : 
I Y  ’  Y  * : 
IZ'Z’ : 

IX*  Y*  : 
IX'Z' : 
I Y*  Z  ’ : 


0.791766E-02 

0.620526E-02 

0.312264E-02 

0. 

-.657281E-03 

0. 


Figure  B-2.  Example  problem  data 
fi  le. 


Figure  B-3.  Example  problem  solution. 
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